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Abstract 

A detailed study of an inhomogeneous dust cosmology contained in a 7-law 
family of perfect-fluid metrics recently presented by Mars and Senovilla is per- 
formed. The metric is shown to be the most general orthogonally transitive, 
Abelian, G2 on S2 solution admitting an additional homothety such that the 
self-similar group H3 is of Bianchi type VI and the fluid flow is tangent to its or- 
bits. The analogous cases with Bianchi types I, II, III, V, VIII and IX are shown 
to be impossible thus making this metric privileged from a mathematical view- 
point. The differential equations determining the metric are partially integrated 
and the line-element is given up to a first order differential equation of Abel 
type of first kind and two quadratures. The solutions are qualitatively analyzed 
by investigating the corresponding autonomous dynamical system. The space- 
time is regular everywhere except for the big bang and the metric is complete 
both into the future and in all spatial directions. The energy-density is positive, 
bounded from above at any instant of time and with an spatial profile (in the 
direction of inhomogeneity) which is oscillating with a rapidly decreasing am- 
plitude. The generic asymptotic behaviour at spatial infinity is a homogeneous 
plane wave. Well-known dynamical system results indicate that this metric is 
very likely to describe the asymptotic behaviour in time of a much more general 
class of inhomogeneous G2 dust cosmologies. 
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1 Introduction 



Solutions of Einstein's field equations for dust with the energy-momentum tensor 

T a(i = pU a Uf3, (1) 

where p and u are the energy-density and the velocity vector of the fluid, respectively, 
are adequate for describing cosmological models at late times, i.e. when the electro- 
magnetic radiation becomes dynamically negligible and the evolution of the universe is 
dominated by matter. The standard model describing the geometry of the universe at 
late times is the spatially homogeneous and isotropic Friedman-Lemaitre-Robertson- 
Walker metric with vanishing pressure. However, since the universe is not exactly 
isotropic and homogeneous, it becomes necessary to consider also less symmetric dust 
cosmologies in order to describe those irregularities. The first natural generalization are 
the spatially homogeneous Bianchi (including locally rotationally symmetric - LRS-) 
cosmologies |l] (see also @ for more recent developments). 

The next step forward is generalizing the models so as to describe inhomogeneous 
cosmologies. The simplest case here is the class of inhomogeneous models admitting 
a three-dimensional group of isometries acting on spacelike two-surfaces (G3 on S2 
models). That class of dust solutions includes the spherically symmetric case (Tolman- 
Bondi metric) first found by Lemaitre || and Tolman Q (see also Krasihki's review 
book || for a comprehensive discussion on this topic) and the plane and hyperbolic 
symmetric ones first discovered by Ellis ||. The study of more general inhomogeneous 
models has been much less systematic and only very partial and incomplete results 
are known. Among them, we find the remarkable Szekeres JF| class of inhomogeneous 
irrotational dust solutions with Petrov type D and vanishing magnetic part of the 
Weyl tensor (along the fluid congruence). A detailed analysis of the physical and 
mathematical aspects of this class of solutions can be found in H and references therein. 

In a systematic approach for investigating inhomogeneous dust cosmologies, the 
next step is the study of spacetimes admitting two linearly independent spacelike Killing 
vector fields, the so-called G2 cosmologies. Very few solutions are known within this 
class. The first one is contained in a 7-law family of solutions obtained recently in 
|8[]. The solution was given up to a coupled system of non-linear, ordinary differential 
equations. As we shall see below, this metric has interesting mathematical and physical 
properties and will be carefully studied in this paper. Very recently, the first completely 
explicit G2 dust solutions have been obtained ||. It is worth mentioning here that the 
situation is completely different when the spacetime is stationary and axially symmetric 
(so that the isometry group is still two-dimensional and Abelian, but acting on timelike 



surfaces). See [10] for a review of dust solutions for this case. 

In this paper we study in detail the inhomogeneous dust solution of Einstein's field 
equations contained in the 7-law family of perfect fluids given in ||. Its maximal 
isometry group is two-dimensional, Abelian and acting orthogonally transitively on 
spacelike orbits. The metric is non-diagonal and turns out to be self-similar, i.e. the 



spacetime possesses also a homothetic Killing vector (see [ I I] for a definition). The 
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self-similar group H3 (spanned by the Killing vector fields and the homothetic Killing 
vector) acts on timelike hypersurfaces and the fluid velocity is tangent to its orbits (we 
will use the term tangent self-similar solution to denote those solutions for which the 
fluid velocity is tangent to the self-similar orbits). As we shall see below, this metric is 
the most general orthogonally transitive G2 tangent self-similar solution such that the 
H 3 is of Bianchi type VI. Furthermore, no analogous solutions for Bianchi types I, II, 
III or V exist. Since Bianchi types VIII and IX are impossible (they do not contain an 
Abelian two-dimensional subgroup), the only remaining cases which could in principle 
admit self-similar tangent solutions are Bianchi types IV or VII (these cases are not 
considered in this paper). As explained below, the tangent self-similar inhomogeneous 
models are known to be suitable candidates to describe the behaviour of more general 
cosmologies at late times or near the big bang singularity. Therefore, the solution 
studied here is very likely to describe the asymptotic behaviour in time of a large class 
of general inhomogeneous G2 dust cosmologies. The metric is algebraically general, 
sharing this property with the explicit solutions recently found in H (both Szekeres 
and G 3 on S 2 are of Petrov type D). The system of coupled, ordinary differential 
equations is partially solved so that the metric can be given up to one first order 
ordinary differential equation of Abel type and two quadratures. The most relevant 
properties of the solution are obtained by performing a detailed qualitative study of 
the dynamical system associated with the Abel equation. In particular, we show that 
the metric is free of singularities except for the big bang, the spacetime is complete 
both into the future and all spatial directions and the energy-density is positive and 
bounded from above at each instant of time. The number of essential parameters is 
two. 

The plan of the paper is as follows. In section 2, the line-element given in || and the 
differential equations to be solved are written down. The results stating the uniqueness 
of the metric are presented and a summary of known results about the relevance of self- 
similar cosmologies as asymptotic states for more general models is included. We then 
rewrite the metric so that only a first order differential equation and two quadratures 
remain to be solved and we give some properties of the spacetime, like the kinematical 
quantities of the fluid, the energy-density, the Petrov type and the Killing generators 
of the self-similar group, as well as the Bianchi type of the Lie algebra they generate. 
In section 3, a detailed qualitative analysis of the two-dimensional dynamical system, 
including a description of the topology of the phase space, the set of equilibrium points 
and its character and a proof of the non-existence of periodic orbits is performed. 
We then consider the asymptotic behaviour of the solution at spatial infinity which 
turns out to be a homogeneous vacuum plane wave of Petrov type N. The coordinate 
transformation which brings this plane wave metric into the standard form is then 
given, thus allowing us to partially interpret the original coordinates in which the 
dust metric is written. Finally, the limiting cases for the family are analyzed. In 
section 4 we describe how this dust metric fits into the general framework developed 
by Hewitt and Wainwright to describe Abelian G2 perfect-fluids with a 7-law equation 
of state. This is particularly interesting due to the special character of this solution 
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as an equilibrium point of the dynamical system in terms of expansion normalized, 
dimensionless variables. 



2 General dust solution 

Finding perfect-fluid solutions of Einstein's field equations for G<i Abelian cosmologies 
is a difficult task and several simplifications have been used in order to obtain par- 
ticular solutions. One of the most fruitful ansatze has been assuming that the metric 
coefficients in comoving coordinates are products of functions of only one variable (the 
so-called separable cosmologies). Within the class A(ii) in Wainwright's classification 



12| for G2 Abelian cosmologies, the separable case with a 7-law equation of state (i.e. 



p = 7p where p is the pressure of the fluid) has been studied in [pi]. From results in 
that paper it is easy to infer that the particular case 7 = implies vacuum and so no 



dust solutions is possible. The separable diagonal case was studied in |L4| (see also [IS 



for the stiff fluid solutions and the pioneering work in |T(| where the solutions with the 



three-spaces orthogonal to the fluid velocity being conformally flat were found). Again, 
no dust solutions are contained in this family. The next case in complexity (class B(i) 
in Wainwright's classification) has been recently analyzed in || where the complete list 
of solutions has been given (for some families a set of ordinary differential equations 
remains to be solved). One of these families describes a perfect fluid with 7-law equa- 
tion of state. The particular case 7 = (dust) is now possible. From the results in || 
it is not difficult to see that the most general dust solution for spacetimes admitting 
an Abelian maximal two-dimensional isometry group acting orthogonally transitively 
on spacelike surfaces with separable metric coefficients in comoving coordinates can be 
written in the form 

ds > = e 2 b i (_ d p + d ^ + GPe o+i)t dy 2 + ^ e (b-Di ^ dz + We t dy y } (2 ) 

where b is an arbitrary non- vanishing constant and G, P and W are functions of x which 
must satisfy the ordinary system of differential equations (prime denotes derivative with 
respect to x) 

G' P' WW n 
P' 2 W' 2 w 2 

^ + ^-7^ = o> (3) 

W ' 2 W2 ,2 1 n 

" pT" pT + 6 2 "1 = 0. 

The symmetry assumptions in deriving this metric involved only a maximal isometry 
group of dimension two. However, it turns out that the metric possesses also a 
homothetic Killing vector (in fact, this also holds for the full class of solutions with 
7-law equation of state in PI). The three-dimensional homothety group H3 belongs 
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to the Bianchi type VI and its orbits are timelike everywhere with the the fluid flow 
tangent to them. 

Spacetimes admitting a three-dimensional homothety group have been investigated 
by several authors in the last few years (see e.g. |HJ and references therein). In 
particular, a recent paper by Carot et al [18|] gives the explicit time dependence for 
this kind of metrics. Furthermore, in |T7j it is stated that tangent self-similar perfect 
fluids of Bianchi types I and II have to satisfy p + 3p = thus excluding dust solutions. 
Combining their results with the uniqueness of the metric (Q) as the only dust solution 
with separable coefficients in comoving coordinates, the following proposition holds 

The line- element ([Jj satisfying the equations (Qj is the most general dust solution 
of Einstein's field equation for spacetimes of class B in Wainwright's classification of 
inhomogeneous cosmologies, admitting one homothety and such that the fluid flow is 
tangent to the self-similar group orbits and the Bianchi type of the three-dimensional 
homothetic group is VI. Furthermore, the analogous cases in which the self-similar 
group is of Bianchi types I, II, III and V are impossible. 

Additionally, Gi Abelian cosmologies with an additional homothetic vector field 
such that the three-dimensional self-similar group H 3 is of Bianchi types VIII and 
IX are clearly impossible because these Bianchi types do not contain Abelian two- 
dimensional subalgebras. In principle, G2 dust solutions with an H 3 of Bianchi types 
IV and VII could admit tangent solutions but none of these algebras admit diagonal 
metrics (unless the spacetimes possess more than two isometries). Furthermore, the 
diagonal limit of the metric (|2|) (obtained for W = 0) can be trivially seen to be either 
Minkowski or a Bianchi cosmology. Thus, the following result also holds 

No dust solutions of Einstein's field equations exist for diagonal inhomogeneous 
cosmologies (B(ii) in Wainwright's classification) admitting a homothetic symmetry 
such that the fluid flow is tangent to the self-similar orbits. 

Tangent self-similar cosmologies play a prominent role in describing the asymptotic 
behaviour (either near the big bang singularity or at t — > 00) for more general cosmo- 
logical models. This result is one of the main consequences of the dynamical system 
analysis of cosmological models which has been undertaken in recent years (see f|| and 
references therein). The dynamical system approach consists of writing Einstein's field 
equations as first order partial differential equations for some kinematical quantities 
associated with an orthonormal tetrad which must be chosen appropriately (T^|. Then, 
an expansion normalized set of variables is introduced and a qualitative analysis of the 
resulting dynamical system is performed. This approach for analyzing cosmological 
models started by considering the simplest cases, namely homogeneous cosmologies in 
which the Einstein field equations become ordinary differential equations and there- 
fore the associated phase space is finite dimensional. The mathematical theory for 
describing finite dimensional dynamical systems is well-developed and a good number 
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of properties of the solutions can be derived without the need of solving the equations 
explicitly. One of the most relevant results which were proven in that context f20fl , 



2T| is that all equilibrium points of the dynamical system associated with non-tilted 
Bianchi cosmologies (i.e. such that the fluid velocity is orthogonal to the group or- 
bits) satisfying a 7-law equation of state correspond to self-similar cosmologies, i.e. 
spatially homogeneous spacetimes admitting a homothetic Killing vector. Thus, this 
class of self-similar Bianchi models plays a relevant role as asymptotic states for more 
general cosmologies. 

For the case of G2 cosmologies, the phase space becomes, in general, infinite dimen- 
sional since the tetrad form of Einstein's equations are first order partial differential 
equations. The mathematical theory in this case is less developed and the results 
which can be obtained are more restricted. However, Hewitt and Wainwright p3] re- 



formulated the field equations as a dynamical system for orthogonally transitive G 2 



perfect fluid cosmologies with a 7-law equation of state. In ]22| it has been proven 
that the equilibrium points of that dynamical system also correspond to self-similar 
models. They admit a three-dimensional self-similar group H s (spanned by the two 
Killing vectors and one homothetic Killing vector) acting on timelike hypersurfaces. 
The fluid velocity vector is tangent to the three-dimensional similarity orbits. These 
two results suggest that tangent self-similar models are important in describing the 
asymptotic behaviour of a large class of inhomogeneous cosmologies. Thus, the unique 
(for the Bianchi types I, II, III, V and VI of the H3 group) tangent self-similar solution 
(H) is very likely to be relevant for describing the asymptotic state of general dust 
cosmologies. 

All these results above show that the metric ([|) is privileged from a mathematical 
point of view, and, thus, it deserves a more detailed analysis. Furthermore, as we shall 
see below, it turns out that the metric is also well-behaved from a physical viewpoint. 
Let us then proceed with the analysis of the metric (H) and field equations @. First 
of all, it is immediate to see that the transformation 

P W 
« — > —b, t < — > —t, z < — > y, P < — ► — - — — — -, W 



P 2 + W 2> P2 + W 2 

leaves the metric (0) invariant. Thus, the constant b can be restricted to be positive 
without loss of generality. Furthermore, the energy-density of the fluid is well-behaved 
only for b > 1 (otherwise it is negative everywhere) so that we can restrict the analysis 
to b > 1 and introduce another constant a defined as 

1 

sin a = -, cos a; > 0. 
b 

The system of differential equations we have to deal with constitutes (after substituting 
obtained from the first equation in (|3]) into the other two equations) a non-linear 
coupled system of two ordinary differential equations for two unknowns. One of the 
equations is of second order and linear in the highest derivative and the other is of the 
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first order but quadratic in the highest derivatives. Thus, the system is very compli- 
cated and it must be simplified. To that end we rewrite the metric, after redefining 
the coordinates t and x and rescaling y and 

ds 2 = -dt 2 + t 2 ^- + GPt 1+sina dy 2 + ^-t l - sina (dz + Wt sina dy " 
L(x) P v 

where L(x) is an arbitrary non- vanishing function at our disposal. This function has 
to be chosen appropriately in order to simplify the field equations without introducing 
unnecessary coordinate singularities in the metric. This is not a trivial task and a 
careful investigation of the equations (|3|) is involved. It turn out that the choice 

r/ \ 2 2 2 / 1 / a \ 

L{x) = sin a + cos a cos I — - — I , (4 J 

allows for an explicit integration of the function W and decomposes the remaining 
field equations into a first order ordinary differential equation for the function H and 
two quadratures for G and P respectively, which is a substantial simplification of the 
problem. The metric can be written as 

ds 2 = -dt 2 + — +t 1 + si ™GPdy 2 

sin 2 a + cos 2 a cos 2 ( 

+t 1 " stoa - ( dz + —t sina P cos ( -) dy] . (5) 
P V sin a V 2 / J 

The first order differential equation for the function H is best analyzed when written 
as an autonomous dynamical system in two dimensions. It reads, explicitly, 

dH 

— — = 2cosQ, 
dx 

dQ _ sin Q sin H 

—— = cosQ -\ 5 . (o) 

dx 2 tan a + 1 + cos H 

where we have introduced a new function Q(x). The change Q = arctan [M(H)\ 
transforms this dynamical system back into the first order ordinary differential equation 



dM 1 



dH 



(l + M : 



1 + 



M sin H 



2 tan 2 a + 1 + cos H 



which is an Abel equation of the first kind f23|| . It does not belong to the known inte- 
grable cases f24|. Although apparently simple, solving this Abel equation is still very 



difficult and no particular solutions have been found. However, a qualitative analysis 
of the dynamical system (0) is enough for a good understanding of the solutions. This 
will be performed in the next section. 
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The other two metric coefficients can then be calculated by the quadratures 



= sinQ, (7) 

P x 2 sin a sin Q + cos 2 a cos Q sin H 
P 2 sin 2 a + cos 2 a(l + cosH) 

The number of essential parameters in the family ([5]) is two; the explicit constant a in 
the line-element and the initial condition in the dynamical system (^) . The integration 
constants in G and P are superfluous and can be set equal to one by trivial redefinitions 
of y and z. 

Let us now describe some of the properties of the dust solution we are considering. 
The fluid velocity vector is u = and the energy-density reads 



cos 2 a sin 2 



P= ^ (9) 

which is positive everywhere and bounded from above at any instant of time t. The 
metric @ has a big bang singularity at t — 0, where the energy-density blows up. 
Regarding the kinematical quantities, the fluid velocity is obviously geodesic and irro- 
tational and the expansion is spatially homogeneous and reads 

•4 

showing that the fluid never recollapses. The shear tensor is highly anisotropic and the 
shear scalar is 

a 4-3cos 2 asin 2 (2) 
Vafi* ~ g£j • 

All fluid elements are expanding and a comoving observer would measure red shift in 
every direction. Regarding the Petrov type, the metric is algebraically general (Petrov 
type I) except for the particular case a — | when it degenerates to conformally flat. 
In this case the energy-density also vanishes and the metric becomes flat. The 
magnetic part of the Weyl tensor along the fluid velocity vector in the spacetime (|5]) 
is non-vanishing. Except for some particular subcases described below, the metric (|5]) 
has the only two Killing vectors £2 = d y and £3 = d z . The homothetic Killing vector is 

-> d (1 — sin a) d (1 + sin a) d , . 

and the three-dimensional homothetic Lie algebra spanned by these three vectors is of 
Bianchi type VI^, where h is given by 

h = \- = -b 2 . 

sin a 



This dust solution can be matched to a vacuum solution with two commuting 
spacelike Killing vectors (which can therefore be interpreted as a gravitational wave), 
but the vacuum metric cannot possess an additional homothetic Killing vector any 
longer. Indeed, it is easy to see that the matching between the dust solution and any 
vacuum metric must be performed across one of the homothetic group orbits in the 
dust solution. If we assume that the homothetic Killing vector extends to the vacuum 
region, the time dependence of the vacuum metric becomes determined everywhere and 
must coincide (due to the continuity of the first fundamental forms on the matching 
hypersurface) with the time dependence in the dust region. Thus, the vacuum metric 
must be a homogeneous plane wave. However, the remaining set of matching conditions 
forbids the matching of a homogeneous plane wave and the self-similar dust cosmology 
at finite distances. The field equations for the vacuum solution matching the metric 
(H) are still quite complicated, and no explicit solution has been found. 

3 Phase space analysis and limiting solutions 

In this section we will analyze the spatial behaviour of the dust metrics (|5|). To this 
aim, we have to study the dynamical system (|B|). This is a two-dimensional dynam- 
ical system, depending on an arbitrary parametric constant a, which possesses the 
two discrete symmetries Q —>■ Q + 2n and H — ► H + 2n. There is another dis- 
crete symmetry {Q —>■ —Q,H — > —H}. Thus, it is sufficient to examine the domain 
U = (— 7r, 7r] x (— 7r, 7r] in the phase space diagram for H and Q and use the above 
periodicity. More strictly, the phase space of the dynamical system (|D is topologically 
a two-dimensional torus S 1 x S 1 which can be described by the coordinates (H, Q) EU 
by identifying the two vertical and the two horizontal boundaries respectively. The 
qualitative behaviour of the dynamical system is independent of the value of the pa- 
rameter a and no bifurcations occur. A plot for the phase space portrait with a typical 
value a = it/A (which will be assumed for all Figures) is given in Fig.l. In order to 
understand the behaviour of the solutions we must consider the set of fixed points of 
the dynamical system. There are four fixed points; two of them are saddle points, one 
is an attractor and the remaining one is a repellor. They are given by 



The two solutions corresponding to the two-saddle points are equivalent and represent 
a spatially homogeneous self-similar cosmology (with the three-dimensional isometry 
group of Bianchi type VI). The two solutions corresponding to the attractor and repellor 
are also equivalent and represent vacuum spacetimes. Since the homotopy group of the 
torus is non-trivial (it is isomorphic to the additive group Z), any solution starting 
at the repellor and finishing at the attractor can be classified by the number of times 
the trajectory goes around the torus before reaching its endpoint (i.e. they can be 




(11) 



(12) 
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classified by the homotopy class they belong to). In order to understand the possible 
behaviour of the solutions, it is convenient to draw the set of trajectories finishing 
or starting at the saddle points. The plot is given in Fig.2. These curves divide the 
phase space diagram into basins of attraction which give a clear description of the 
behaviour of the solutions. In particular, we can conclude from this diagram that 
no periodic solutions of (^j) exist. This follows from the fundamental fact that two 
different orbits in any dynamical system cannot intersect anywhere and that any closed 
orbit which is contractible to a point must contain an equilibrium point in its interior. 
Thus, all the solutions of (|]), except those starting and finishing on the saddle points, 
approach asymptotically the repellor (for x — > — oo) and the attractor (for x — > +00). 
Regarding the different homotopy classes the solutions can belong to, it follows from 
Fig 2. that two different cases are possible. Either the solutions go from the repellor to 
the attractor without going around the torus, or they go around it at most once. In the 
second case, the solutions can belong to two different basins of attraction, which can 
be distinguished by the number of times the boundary of U is crossed. This boundary 
can be crossed either once (in the H direction) or twice (one in the H and one in the 
Q direction). A typical solution for each one of the three possible cases is shown in 
Fig. 3, together with the energy- density p, which oscillates with rapidly decreasing 
amplitude. 

A numerical analysis of the solutions of @ shows that the asymptotic approach 
to the stable fixed points ( |1"2"D is quite fast and that the solutions oscillate around 
those points in their asymptotic approach. Since every crossing of the value H = 
corresponds to a vanishing value for the energy- density (see (|9|)), it turns out that the 
energy-density profile of the solutions is oscillating in the positive range, reaching the 
zero value at a discrete number of points and with a rapidly decreasing amplitude (see 
the energy-density plots in Fig. 3 for examples of this behaviour) 

Let us now study in detail the line-elements corresponding to the fixed points, since 
this will allow us for a partial interpretation of the coordinate system {£, x, y, z} in (|5]). 

The asymptotic vacuum solution is represented in the coordinate system (||) by the 
line-element 

ex\l+sina » 2 , /_t sina ( j . COS Q, 



ds 1 = -di 2 + t 2 dx 2 + (te ex Y +sma dy 2 + (te tx y- sma dz + (te ex ) sma dy (13) 

V sin a / 

where e = +1 (e = —1) for x — > +oo (x — > — oo). This vacuum spacetime is a 
homogeneous plane wave with a six- dimensional isometry group acting transitively on 
V4. The Petrov type is N and the only repeated principal null direction is I — td t — ed x . 
Thus, the inhomogeneous dust metric contains a matter and a radiative part. The 
radiative part dominates at large distances in the spatial direction perpendicular to the 
two-planes spanned by the isometries. In a cosmological context it could be interpreted 
as a gravitational wave background superposed to the expanding dust. This asymptotic 
behaviour at spacelike infinity seems to be quite common in this kind of situations for 
perfect fluids. Hewitt et al [|25[] performed a dynamical system analysis of tangent 



self-similar diagonal G2 perfect fluids with a 7-law equation of state. They concluded 
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that the spacelike asymptotic behaviour can be either matter dominated (i.e. with a 
non-zero value for the energy- density at spacelike infinity) or vacuum dominated (in 
which the energy- density tends to zero at infinity). In this second case, the vacuum 
asymptotic behaviour is also a homogeneous plane wave. 

The very neat behaviour of the solution at spatial infinity can provide us with 
a partial interpretation for the coordinate system {t,x,y,z}. In order to do so, we 
consider the coordinate change which brings (^) into the standard form for plane 
waves [|TJ 



ds = -2dudv + 2 



A(u) [Z 2 -Y 2 ) + 2B(u)YZ du 2 + dY 2 + dZ 2 . (14) 



It can be easily seen that this coordinate transformation is given by 

— l + sin ol 

y = _ M - i± f iii [Y sm((3 -a) + Z cos(/3 -a)], z = - — - — (Y sin (3 + Z cos (3) , (15) 

sin a 

te £X = u, te~ ex = 2v- — 

2u 

where 



Z 2 - Y 2 ) sin(2/3 - a) - 2YZ cos(2/? - a) + Y 2 + Z 2 



COS Ol , . 

= — logu. (16) 

The two function A{u) and B{u) in ( |I^) read 

. . . cos a . s cos a . . 
A(u) = -cos(2p-a), B(u) = r sin2/J-a 

So, near spatial infinity a partial interpretation of the coordinate system (|5[) is possible. 
The coordinate x measures a spacelike distance along the direction of propagation of 
the plane wave. The coordinates y and z are straight lines lying in a plane tangent to 
the null hypersurface corresponding to the plane wave front. They are not orthogonal 
lines, the angle between them being directly related to the parameter a. The lines 
y = const, and z = const, are scaled by different units, so there is no rotational 
symmetry in the corresponding plane. Moreover, the straight lines associated with y 
and z are rotated by an angle f3 (which depends on u and therefore is different for 
different wave fronts) with respect to the orthogonal straight lines Y, Z associated 
with the plane of symmetry of the wave. It follows then that the coordinates y, z are 
rotating with respect to the homogeneous plane wave background. 

The solution corresponding to the saddle points is a diagonal spatially homogeneous 
spacetime of Bianchi type VI^_ sin2Q ,^. The line-element ([5]) takes the form 

ds 2 = -dt 2 + t 2 dx 2 + {te x ) 1+sina dy 2 + {te- x Y~ sina dz 2 (17) 

which possesses the additional killing vector field 

d 1 + sin ad 1 — sin a 
^ = d-x-^2— y d-y + ^2— zdz 
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and the energy density is p = co % - . This metric was first obtained by Collins p6[ and 



appears listed in the catalogue of self-similar solutions by Hsu and Wainwright [27 



This solution also appears as a saddle point in a dynamical system investigation of 



tangent self-similar, G2 diagonal perfect fluids with a 7- law equation of state | 25fl . 

For a = I the line-elements ([TJ]) and (|I7J) are exactly the same. Moreover, the 
expressions for A(u) and B(u) vanish when a = ~ and therefore the metric becomes 
flat. The transformation ( |l5l) particularized to a — f simply reads 

Y Y 2 

y = — , z = Z, te x = u, te x = 2v , 

u u 

which can be rewritten using the standard orthogonal coordinates T, X, Y, Z for Min- 
kowski spacetime as 

t = J T 2_ X 2_ Y 2 e = (T + X) V2Y 

The four velocity u = — dt is transformed to 

u = 1 = (TdT + FdY + XdX) . 

VT 2 - X 2 - Y 2 V ; 

Let us finish the discussion on the solution (||) by considering another limiting case. 
The family (||) is defined for every value of a except when sin a = (this corresponds 
to b = 00 in the metric (|2])). It is therefore convenient to consider the limiting case 
sin a — > 0. It can be easily seen that rescaling the coordinates in order to reabsorb 
the diverging behaviour in the coefficient dydz in (0) leads necessarily to a degenerate 
metric. Therefore, the only possibility is that cos (-rf) tends to zero when sin a — > 0. 
Since the metric coefficient in dx 2 then diverges, it is necessary to perform a coordinate 
transformation to obtain a regular metric. Thus, we define 

cos ( — ^ = 5(r) sin a, . = dr 

which brings both the metric and the field equations into a form allowing the limit 
sin a — > 0. The field equations turn out to be completely integrable and the resulting 
line-element can be written as 

ds 2 = -dt 2 + t 2 dr 2 H — dy 2 + tcoshr ( dz + tanhr dy] , (18) 

cosh r \ J 



which can be transformed into the particular case sin a = in ( [LTD by the coordinate 
transformation 

y+z y- z 

r — > x, — = > y, — = > z. 

V2 V2 
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Thus the solution ( |TTD is the correct spatially homogeneous limiting case for all values 
of a. 

As we have emphasized above, the dust model we study in this paper corresponds 
to an equilibrium point in the dynamical system approach for G 2 Abelian perfect-fluids 
developed by Hewitt and Wainwright |[22]| . It is, therefore, interesting to analyze how 
this solution fits into that framework. This will be the task of the next section. 



Description of the solution using the tetrad ap- 
proach 



Throughout this section, we shall adopt the same notation and definitions as in |22[ 



One of the main results in that paper is that every equilibrium point of the dynamical 
system they construct (in terms of expansion normalized, dimensionless connection 
coefficients) corresponds to an orthogonally transitive G 2 perfect-fluid (assuming a 
linear equation of state) which admits a homothetic vector field lying in the three- 
plane spanned by u and the two Killing vectors (i.e. a tangent self-similar solution). 
Since these conditions are met by the metric we analyze, we can assume OqX = where 
X represents any dimensionless variable. Furthermore, the fluid is dust and hence the 



four velocity is geodesic. Imposing these conditions on the dynamical system in |22l , 
we readily obtain 

,4 ^ = 4, -0, 

which, in particular, imply that do and d\ commute. Thus, every tangent self-similar 
dust solution has constant and positive deceleration parameter q and the metric is 
never inflating. The defining equations for q and r provide now 

Q = - - 12 (t\ 2 + E 2 3 ) , A =12 (t 22 N 23 - £ 23 A> 22 

indicating a very simple relationship between the density parameter Q and the dimen- 
sionless shear scalar. Using all this information, the dynamical system splits into a set 
of ordinary differential equations 

d 1 R = -18(t 22 N 22 + £ 23 N 23 ), 
djl 23 = -2 (iV + S 22 + RN 23 ) , djl 22 = 2 (iV + £ 23 - RN 22 ) , 

d!iV 2 3 = 24A> 23 (£ 22 iV 2 3 - £ 23 A> 22 ) - 2iV + A> 22 + 2Rt 23 - ^S 22 , (19) 

d!N 22 = 2AN 22 {t 22 N 23 - t 23 N 22 ) + 2iV+iV 23 + 2i?S 22 + ^S 23 , 

and the algebraic constraint 

1 - 16(A> 2 2 2 + A> 2 2 3 + S 2 2 + E 2 3 ) + 16 2 (S 23 iV 22 - £ 22 A> 23 ) 2 = 0, 
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which is a first integral of Ql9]). It is a matter of simple calculation to check that 

R 2 = 9(£ 2 2 + £ 2 3 ) + k (20) 

(where k is an arbitrary constant) is also a first integral of (|H?|). The sign of this constant 
can be related to the Bianchi type of the homothetic algebra. It turns out that k < 
corresponds to Bianchi VI, k = to Bianchi IV and k > to Bianchi VII. Thus, the 
dynamical system is six-dimensional (in the variables R, £22, ^23, ^22, -V23 and N + ) 
with two polynomic constraints. In order to simplify it further we take advantage of 
the rotational freedom in the e2, e 3 plane (the rotation angle (f) must be constant both 
along the group orbits and along the integral lines of the fluid velocity in order to 
comply with all the conditions we already imposed). Under such a transformation, R 
behaves as an scalar, N + transforms as iV + — > N + + di<f>, Sab and Nab as symmetric 
rank two tensors. The Lie algebra we study in this paper is of Bianchi type VI and 
therefore k is strictly negative, which will be assumed from now on. Hence, the first 
integral Q2"DD implies R? < 9(X 2 2 + ^23), which allows us to fix so that the condition 

R = -3S 23 (21) 

is fulfilled. Notice that this choice cannot be imposed in the case of Bianchi type VII 
(for which k > 0). The choice ( pl~D immediately implies E22 — :y y^- The non-negativity 
of the energy-density demands < 1 and we can set 4v ^~^ = sin a. Furthermore, 

the differential equation for E22 above fixes iV + = —3^22- Inserting all this information 
into the dynamical system, we obtain the following set of equations 

3 

di£ 2 3 = 2 ( 4S 23^V23 + sinaiV2 2 ) , 
9^23 = 6 (iY 2 2 3 sina - ^iV^s + N 2 22 - t 2 23 ) - (22) 



diN 22 = 6N 22 [(sin a - 1) N 23 - AN 22 t 23 
and the constraint 



3 (1 - sina) S 23 

T~ _ 



cos 2 a (l - 16iV 2 2 3) - 16S^ 3 (l - 16iV 2 2) - 16iV 2 2 (N 22 + 8V23S23 sin a) = 0. (23) 

Hence, the dynamical system is two-dimensional as the motion happens on the surface 
S defined by (|23|). In order to obtain an adequate description of the problem we should 
find appropriate global coordinates on the surface S such that the dynamical system 
(^) takes a simple form, which is not a trivial problem. Since we have already analyzed 
the solution from the coordinate perspective, it is a matter of simple calculation to 
obtain which coordinates in 5* provide the dynamical system studied in section 2. Since 
d and d± commute we can introduce coordinates t and x such that d = 3t/2d t and 



d\ = 3/2JL(x)d x , where L(x) is an arbitrary non-vanishing function. Parametrizing 
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S 2 3, ^23 and N 2 2 as 



cos a cos (f) —2 sin a sin Q — cos 2 a cos Q sin H 

£23 = ^ J ^23 — j ; 

'sin 2 a + cos 2 a cos 2 ( yj 



22 



cos a 
4 sin a 



I cos ( — ^ A^23 + cos Q sin y sin 2 a + cos2 a cos2 f^" 



in terms of H an Q (which are global coordinates on the surface S) the dynamical 
system ( p2|) takes exactly the form @ whenever we also choose the free function 
L(x) as in (|4]). Clearly, finding this parametrization of S by simply inspecting the 
dynamical system (|22|) is very difficult. If we had started from the tetrad approach we 
would probably have found a different parametrization of S but the description of the 
problem would have been essentially equivalent. Actually, it is remarkable that both 
methods lead quite naturally to a two-dimensional dynamical system. This implies that 
both the level of final understanding and the amount of the required work are similar for 
the two approaches. Performing the comparison of them is, however, very interesting 
because the dimensionless variables of the tetrad formulation have a very clear physical 
interpretation. Furthermore, the tetrad approach is coordinate independent and hence 
it characterizes the solutions intrinsically. Finally, a combined use of both methods can 
provide us with a powerful tool for solving more difficult problems. In particular, the 
analysis of Bianchi type VII tangent self-similar dust cosmologies has proven to be very 
difficult following the coordinate approach. It is possible that using and combining the 
information from the coordinate and the tetrad methods can help solving this case. This 
matter is now under current investigation and the results will be reported elsewhere. 
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Figure 1: Phase space portrait of the dynamical system @ for a = W. The fixed 
points are marked by circles 




Figure 2: Phase space diagram for the dynamical system (J|) showing the set of trajec- 
tories starting or finishing on the saddle points Ok, ~) and (jr, (which correspond 
to four points in the diagram due to the identification of the two vertical lines). The 
direction of increasing x for each solution is shown. 




Figure 3: Phase space portraits and energy-densities for three different solutions of 
the dynamical system with a = tt/A. The thin-line solution corresponds to the initial 
values (5(0) = | and H(0) = 0. It never crosses the boundary of U and it does not wrap 
the torus. The medium-line solution corresponds to the initial values Q(0) = and 
H(0) = f 7T. This solution wraps the torus once and the boundary of U is crossed once 
in the H direction. The thick-line solution corresponds to Q(0) = |vr and H(0) = |7r. 
It wraps the torus once and the boundary of U is crossed twice, once in the H and 
once in the Q direction. 



